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Abstract 

Head-to-head and tail-to-tail 180° domain-walls in a finite isolated ferroelectric sample are theo- 
retically studied using Landau theory. The full set of equations, suitable for numerical calculations 
is developed. The explicit expressions for the polarization profile across the walls arc derived for 
several limiting cases and wall- widths are estimated. It is shown analytically that different regimes 
of screening and different dependences for width of charged domain walls on the temperature and 
parameters of the system are possible, depending on spontaneous polarization and concentration 
of carriers in the material. It is shown that the half-width of charged domain walls in typical per- 
ovskites is about the nonlinear Thomas-Fermi screening-length and about one order of magnitude 
larger than the half- width of neutral domain- walls. The formation energies of head-to-head walls 
under different regimes of screening are obtained, neglecting the poling ability of the surface. In 
the nonlinear regimes of screening, this energy is equal to the energy necessary for the creation of 
electron-hole pairs in the amount sufficient to screen the spontaneous polarization, which is pro- 
portional to the band gap of the ferroelectric. It is shown that either head-to-head or tail-to-tail 
configuration can be energetically favorable in comparison with the monodomain state of the fer- 
roelectric if the poling ability of the surface is large enough. If this is not the case, the existence of 
charged domain walls in bulk ferroelectrics is merely a result of the domain-growth kinetics. For- 
mation energies of the other possible states: multidomain state with antiparallcl domains separated 
by neutral walls and the state with the zero polarization were compared with the formation energy 
of the charged domain wall. It was shown that, at large enough sample thicknesses, a charged do- 
main wall can be energetically favorable in comparison with the states mentioned above. This size 
effect could explain why charged domain walls were observed experimentally in bulk lead titanate 
but not in barium titanate. The results obtained for the case of an isolated ferroelectric sample 
were compared with the results for an electroded sample. It was shown that charged domain wall 
in electroded sample can be either metastable or stable, depends on the work function difference 
between electrodes and ferroelectric and the poling ability of the electrode/ferroelectric interface. 

PACS numbers: 77.80.Dj, 72.20.Jv 
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I. INTRODUCTION 

Ferroelectric materials are widely used today in memories, piezoelectric transducers, py- 
roelectric detectors and thin-film capacitorsii^. The structural and functional properties 
of domain walls can substantially influence poling and polarization switching - two most 
important processes in ferroelectrics and their applications. 

Three qualitatively different configurations are possible for 180-degrees domains (Fig. [1]). 
In the first case, the domain wall is parallel to the direction of the spontaneous polarization 
inside the adjacent domains. In this configuration there is no bound charge on the wall. This 
case is well studied and described in details (see e.g. the book by Strukov and Levanyuk^). 
In the other two configurations, the domain wall is not parallel to the direction of the 
spontaneous polarization. The polarization can be directed either towards the domain wall 
("head-to- head") or from the domain wall ("tail-to-tail") so that positive or negative bound 
charges are present on the domain wall, respectively. 

Neutral "Head-to-head" "Tail-to-tail" 
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FIG. 1. Neutral and two types of charged domain walls, with the polarization normal to the wall 
surface. Orientation of spontaneous polarization with respect to the wall is shown with arrows. 

Although neutral domain walls are much more common, charged domain walls were 
observed in various ferroelectrics, for example in PbTiOs crystals^"-, Pb[Zr2^Tii_a;]03 (PZT) 
ceramics^, BiFeO^^, and PZT thin films^. 

The properties of charged domain walls can be quite different from the properties of 
neutral domain walls. For example Mokry, Tagantsev, and Fousek^" showed that the com- 
pensation of the polarization charge by free carriers reduces the pressure excerted on the 
wall when an external electric field is applied. This leads to a reduction in the mobility of 
the wall and an increase in switching voltage, even when the mobility of the compensating 
free charge is high. This effect was experimentally observed by Balke et. al^-'^ in BiFeOa thin 
film. 
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In a perfect insulator, unscreened bound charge on the charged domain wall typically 
creates a very large electric field, which effectively shifts down the Curie temperature (for 
the wall perpendicular to the direction of polarization in barium titanate this shift is about 
24000 K.) Therefore, except for a very small spontaneous polarization and a very small angle 
between the vector of spontaneous polarization and the domain wall, for charged walls to 
exist, their bound charge should be almost completely screened by free charges. These free 
charges should be taken into consideration in the theory for charged domain walls. It is 
this situation that will be addressed in the paper. If it is not mentioned specifically, we will 
consider "normal" ferroelectrics, with high permittivity 6{ ^ e^, where is the background 
permittivity. The model considered is suitable only for the case of the big angle between 
the domain wall and the vector of spontaneous polarization. 

"Head-to-head" and "tail-to-tail" domain walls in a sample with metallic electrodes were 
considered in a series of works by Ivanchik, Guro, Chenskii et al.— It was argued that such 
domain walls can be stable depending on the work function difference between the metal 
electrode and the ferroelectric. Below we apply the Landau theory to address the problem 
of charged domain walls in an isolated ferroelectric. Some models of "head-to-head" and 
"tail-to-tail" domain walls in an isolated sample were considered in the past, where the 
wall was stabilized by an inhomogeneous distribution of a dopant.— We will consider the 
case of a homogeneous material. The results obtained below for the internal structure and 
characteristic scales of domain walls can also be directly applied to the case of an electroded 
sample as well. We will revisit the problem of domain wall width because several different 
scales were obtained in the past by different authors. For example, according to the classical 
book on ferroelectric-semiconductors by Fridkin-^, domain wall formation in head-to-head 
configuration is controlled by Debye screening with free carriers and the wall thickness should 
be about the Debye screening length. Another scale was found by Krapivin and Chenskii.— 
We also compare the results obtained for the isolated ferroelectric with the results obtained 
in the past for the electroded samples.— >i^ii^ 

The paper develops as follows: First, we formulate the problem in a way suitable for 
numerical calculations (Sect. HI])- In Sect. IIII Al we discuss characteristic scales of the po- 
larization variation and limits of applicability of the continuos theory. In Sect. IIII Bl the 
problem is formulated in a way, suitable for an analytical treatment. We obtain analytical 
solutions and evaluate domain wall widths corresponding to the different regimes of polar- 
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ization screening by free carriers (Sect. IIII CI and IIII Dh The correspondence between the 
regimes of screening and parameters of the ferroelectric is discussed in Sect. IIII E[ 

In Sect. IIVI we evaluate the energies of isolated samples with " head-to- head" and "tail- 
to-tail" domain walls. Different configurations that can be energetically favorable depending 
on the thickness of the ferroelectric, and size-effects are considered in Sect. El An analysis 
is presented in Sect. IVI| where the situation in isolated ferroelectric sample is compared 
with that in an electroded sample. Numerical estimates and comparison of the prediction 
of the theory with experimental data are presented in Sect. IVIII 

II. GENERAL FORMULATION OF THE PROBLEM 

We consider "head-to-head" and "tail-to-tail" domain walls perpendicular to the direction 
of polarization (as shown in Fig. [1]) in an isolated sample. If the ferroelectric material is an 
ideal insulator, a huge depolarizing field will appear, which will suppress ferroelectricity in 
the sample. As a model, we will consider an infinite plate with thickness L. We consider a 
single component system, where all the variables depend only on one x coordinate, where 
the X axis is perpendicular to the plate surface. 

First, we obtain the system of equations for the isolated ferroelectric element using the 
Landau theory, taking into account the screening charges. For simplicity we consider ferro- 
electrics with a second order phase transition. However, most of the results obtained below 
are valid for ferroelectrics with a first order phase transition as well. 

The equation of state for ferroelectrics with a second order phase transition has the form 
(see e.g.—): 



where a < in the ferroelectric state, /3 > 0, and k is the coefficient of the gradient term in 
the free energy, E is the electric field, and P is the ferroelectric part of the polarization. 
The electrical field and potential are connected by the relation: 

The Poisson equation for one-dimensional case has the form: 




(1) 




(3) 
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where p is the free charge density, D is the electrical displacement, defined as: 

D = ehE + AirP, (4) 

is the background permittivity. 

In the model of electron gas with parabolic spectrum the free charge density p depends 
on the potential (f as follows22"^: 

(5) 

A^c and Ny are the effective density of states in conductive and valence band respectively, 
Ey is the top of the valence band and Ec is the bottom of the conduction band, E-p is the 
Fermi level, q is the absolute value of the electron charge, T is the absolute temperature and 
k is the Boltzmann constant, F1/2 is the Dirac-Fermi integral, defined as: 



00 



y/n J 1 + exp(^ - a) ' 



= — / , , (6) 





Za and Zd are the acceptor and donor valencies, and ta are the fraction of ionized donors 
and acceptors respectively given by: 

where E^ and E^ are the donor and acceptor level respectively, gd is the ground state 
degeneracy of the donor impurity level and g^, is the degeneracy of acceptor level. As we 
will show later, the carriers concentration in a charged domain wall is much smaller than 
the half-band filling one, and so that the use of the parabolic approximation is justified. 

The isolated ferroelectric element is electrically neutral. There is no electric field outside 
the sample, thus the boundary conditions can be formulated as follows: 

^U=-L/2 = /^U=L/2 = (9) 

The other pair of boundary conditions given as (see e.g.—): 

(10) 



x=±L/2 



where C and t] are coefficients in the surface free energy expansion with respect to P. 

Equations (II])- (E]) constitute the full set of equations for the problem. This set of equations 
with boundary conditions Eqs. ([9]), (fTOj) and additional conditions 

P > 0, if - L/2 < s < 

(11) 

P < 0, if < X < L/2 
P < 0, if - L/2 < x < 

(12) 

P > 0, if < a; < L/2 

can be solved numerically to obtain the polarization distribution in the sample with "head- 
to- head" and "tail-to-tail" domain wall respectively. Due to the symmetry of the problem 
boundary conditions can be also redefined as P(0) = 0, -D(O) = and Eqs. ffTUl) at 
X = L/2. 

We define the domain wall formation energy as: 

rL/2 

W= / [$(x) -$o(x)]da; + $surf (13) 

J -L/2 

where $(x) is the free energy density of a system with the domain wall and ^q{x) = ^Pq + 
jPq + $eg("-o) = jPo + '^cg(^o)5 Pq = a/— $cg('^o) is the electron gas free energy at 
equilibrium concentration of carriers in the homogeneous material. Hereafter we will use 
a shorthand "homogeneous concentration" for this way defined carrier concentration. The 
integration is done over the thickness of the sample. The $surf is the surface energy per unit 
area given as follows2ii2^: 

$surf = CPU=-L/2 - CPU=L/2 + ^PX=-L/2 + \P%=L/2 (14) 

The free energy density of the system consists of the free energy of the lattice, the 
electrostatic energy, and the free kinetic energy density of the electron gas $eg: 

The free energy density of the electron gas can be defined as: 

$eg = ^eg - T5eg (16) 
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where Eeg is the kinetic energy density of the electron gas and Seg is its entropy density. 
The kinetic energy can be found by integration over the valence and conduction bands and 
impurities: 

E,g= J ENy{E)f{E)dE + j ENc{E)f{E)dE + z,NMv)E, + z^N^il -U{^))Ea (17) 

— oo Ec 

where Ny{E) and Nc{E) are the density of states in valence and conduction bands respec- 
tively, f{E) is the Fermi function. For the case of parabolic spectrum this relation can be 
transformed to the form: 

(Ef — En. + og5 \ / Ev — E-p — 009 \ , , , , 

— j -NyFy2 [ i^j. ]+zM.{v)E,-ZdNMv)Ed 

(18) 

where Eg is the bandgap and F3/2 defined as: 



00 



2_ r e^dj 

y/n J 1 + exp(^ - a) ' 



Fy.ia) = — I ^ , 1^^^^ (19) 
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The entropy density for Fermi gas can be found as follows: 

Seg = -kJ2 [fi In + (1 - fi) In (1 - /.)] . (20) 

i 

where summation is done over all possible electronic states, fi is the occupation probability 
for the i-th state. 



III. LIMITS OF APPLICABILITY OF THE CONTINUOS THEORY AND ANA- 
LYTICAL SOLUTIONS IN A LARGE CRYSTAL 

A. Characteristic scales and limits of applicability of the continuos theory 

In the previous section we derived the full set of equations, which can be solved numeri- 
cally to obtain the exact solution of the problem. There exist the number of situations where 
approximate analytical solutions to this set can be obtained. Below we will introduce some 
approximations to within a small parameter eh/£f, where ef = 27r/|a| is the contribution of 
the ferroelectric subsystem to the permittivity of the material, which enable us to get these 
solutions. 
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Using Eq. (jlj), the Poisson equation can be presented as: 



dE dP ^ 

+ 47r— = An p. 21 
ox ox 



After differentiation of Eq. (!2T|) and using Eq. ([2]) one can get: 



dx"^ Sb dip eb 
Equation ( l22l) can be rewritten as: 



^2E + }.E = ^-'^—, (23) 

where 

^° = f^^wwr ^^^^ 

This equation contains two self-consistent characteristic scales for the spatial variation of 
the polarization. The big one corresponds to the case, where the first term in Eq. ( l23i) can 
be neglected, and this equation can be presented as: 



dip dx 



E=^- (25) 



Substituting equation of state Eq. ([T]) into Eq. ( j25l) one can get: 



The characteristic scale 5 can be estimated from Eq. fl2^ as: 



^ oi uf, m ^ + (27) 

2\a\{dp/dp) 



where Tc is the correlation radius, defined as: 



2|a| 



(28) 



Hereafter in the estimates, dp/ dip means typical value of dp/ dp inside the region with 
pronounced variation of polarization (e. g. domain wall). 

As we will see later, for any realistic situation 5 ^ Tc, thus from Eq. fl27|) one can get: 



> 1. (30) 



To check the self-consistency of this scale we should show that when the spatial variations of 
the polarization is controlled by this scale the first term in Eq. f l23|) is small in comparison 
with the second one, i.e. 6 Iq. From Eq. and Eq. (12^ we can estimate the ratio of 
these terms to be about: 

Thus our assumption 5 ^ /q is self-consistent. This result is also supported by more rigorous 
calculations as demonstrated in App. |Al using the exact solutions to Eq. fl26|) . 

As the result of the simplification we have passed from a set of the equations equivalent 
to a fourth order differential equation to one second order equation, Eq. (l26l) . At the same 
time we still have 4 boundary conditions, Eqs. (jH]) and f[TU|) . Thus we have excess of the 
boundary conditions. A problem of such kind is a known in theoretical physics, for example, 
in hydrodynamics of liquid flowing in a tube^^. In that case, the problem was solved by 
introducing the boundary layer, the description of which goes beyond the simplified theory. 
We adapt a similar approach by assuming that there exist a thin layer near the surface, where 
our simplified scheme does not work and where the charge density is small in comparison 
with the dP/dx. In this case Eq. ( pTl) together with Eq. ([9]) leads to the following equation: 

e^E + A-kP = 0. (31) 

Using equation of state Eq. ([1]) one can rewrite it as: 

(9^ P Air 
K—^ + aP + I3P^ + — P = 0. (32) 

Neglecting terms about e\,/ei, the characteristic thickness of the boundary layer can be found 
from Eq. (I32l) as follows: 

We can shift boundary condition D = on a distance / from the surface, and the polarization 
changes on this scale to satisfy the second boundary condition. The condition D = means 
that we neglect the charge in the surface layer, that is exactly the assumption we made 
to obtain scale /, thus the solution with scale / is self-consistent. Such approach can be 
supported by numerical calculations. The numerical solution obtained for the case of linear 
screening with classical gas with exact boundary conditions is shown on Fig. [2l Here it is 
seen a deep minimum, which is corresponds to the border of the surface boundary layer. 
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FIG. 2. Numerical solution with boundary conditions Eqs. ([lO]), ([9]), and P(0) = 0, D{0) = 0. The 
small parameter eb/£f = 0.01, S/rc = 4. The dimensionless surface energy coefficients "qrc/n = 
5000, Crc/iKPo) = 5000. 

For typical perovskites Eq. (I33l) implies / ~ ^/K ^ a, where a is the lattice constant. 
Thus the continuous theory is not applicable inside this surface layer. This means that the 
result of / obtained above can be considered as an indication that I is atomically small. The 
numerical solutions obtained for the full set of equations with full boundary conditions will 
also contain the scale smaller than the lattice constant, thus this solution near the surface 
will be outside of the applicability of the continuos theory. Only the solutions with big scale 
and shifted boundary condition can be used. The consideration of the properties of the 
surface layer requires an ab initio approach. For some materials, like week ferroelectrica^ 
or ferroelectrics of order-disorder type far from the phase transition ^b/^f can readily be of 
the order of one or larger. In this case the inequality f pOj) is not valid anymore, and the 
approximation of the boundary layer becomes not self-consistent. However, at the same 
time, the small scale / becomes macroscopic, of the order of the correlation radius, thus 
the solution can be found by solving full system with exact boundary conditions, without 
leaving the range of applicability of the continuous theory. 

Below we will consider normal ferroelectrics (with ^b/^f ^ 1 )• We will not consider the 
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structure of the surface layer and will use macroscopic boundary condition D = 0. 



B. Problem formulation for analytical solution 

Now that the boundary conditions are defined, we can calculate the polarization distri- 
bution in the sample. Equation ( 126|) can be rewritten as: 

aP + /3P' =iK- idp/d^)-') (34) 

This is an equation with respect to P, and it is convenient to reformulate the boundary 
conditions with respect to P instead of D. From Eq. ( 127|) it is immediately follows 6 > r^- 
With Eqs. ([I]), (H]) one can get: 

D = e^ (aP + PP^ + + A-kP (35) 



In out problem \P\ < \Po\ that leads to \PP^\ < \<yP\, the gradient term |k(9^P/(9x^)| < 
|/t(P/r^)| ~ \aP\, thus in the lowest approximation with respect to e^/ef ^ 1, Eq. (!35|) can 
be rewritten as: 

D = AirP. (36) 
So that boundary condition Eq. (jUj) can be presented as: 

P|x=±L/2 = 0. (37) 

The Poisson equation, Eq. ([3]), in this case can be rewritten as: 

dP 

Let us first consider a "head-to- head" domain wall in a large sample, so that the polar- 
ization inside the domains is close to the spontaneous polarization. Thus, neglecting thin 
surface layer, the polarization distribution can be schematically presented as shown in Fig. 
[3l We define the electrical potential as equal to zero in some point inside the domains 
(region 2, Fig. [3]), where the electron concentration is equal to the homogeneous one. For 
the "head-to- head" domain wall (region 3, Fig. [31(a)), for the situation addressed, the prob- 
lem can be formulated as the search of a solution to the set of equations ([I])-® with the 
boundary conditions 

P\-oo = Pq^ P\oo = —Pq- (39) 
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FIG. 3. Schematic polarization distribution in the case of " head-to- head" (a) and "tail-to-tail" 
(b) domain walls. Three different regions are marked for each configuration: 1- surface adjacent 
region, 2 - internal domain region, 3 - domain wall. 



The polarization distribution in the surface adjacent region (region 1, FigJHt^a)) can be 
found from the same equations with the following boundary conditions: 

P|_oo = -Po,P|l/2 = 0, (40) 

P|oo = Po,P|-L/2 = 0. (41) 

This problem is equivalent to the problem for the "tail-to-tail" domain wall (region 3, 
Fig. [3t^b)), which can be found as the solution to Eq. fH5|) with boundary conditions: 

P\-oo = -P^.P\oo = P^. (42) 

Specifically, the polarization profile near the surface of the sample with " head-to-head" wall 
will be a half of the polarization profile in the "tail-to-tail" domain wall. In turn, the solution 
for the "tail-to-tail" domain wall is equivalent to the solutions for the "head-to-head" wall, 
but the screening is provided by holes instead of electrons. To get the results to the "tail- 
to-tail" domain wall, the effective mass and homogeneous concentration of the electrons in 
the results for the "head-to- head" wall should be replaced by the mass and homogeneous 
concentration of holes and vice versa. Further in this section we will consider the specific 
case of a "head-to-head" domain wall. 
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We will consider the limiting cases, when analytical solutions are possible. For the case 
K ^ {dp/d{p)~^ one can rewrite Eq. flM|) as: 

aP + pP' = K^. (43) 

This is the same equation as well known equation for the neutral domain wall. The coefficient 
K, is different for the cases when the polarization gradient is parallel and perpendicular 
to the direction of polarization, but typically this coefficient are of the same order. The 
corresponding solution for the "head-to-head" wall is: 

P = -Potanh(|-). (44) 

In the opposite case k ^ {dp/dip)~^, i.e. 6 Vc this equation can be simplified as: 

|(a.P + /,p3)^-|^. (45) 

We will address four qualitatively different limit situations, where analytical dependence 
p{(p) can be obtained and used in getting the explicit form of dp/ dip. First of all, depending 
on the spontaneous polarization and the homogeneous concentration of electrons, the elec- 
tron gas in the conduction band inside the domain wall can be degenerate or non-degenerate. 
In addition, linear and nonlinear regimes of screening are possible. We should mention that 
the gas inside the wall can be degenerate, even if the electron gas in the conduction band 
in homogeneous ferroelectric is non-degenerate. We will consider the cases of intrinsic semi- 
conductor or fully ionized doping impurities. 

C. Screening with classical electron gas 

In the case of non-degenerate electron gas which obeys classical statistics, Eq. (jS]) can be 
rewritten in the form: 

-Ueoei-T + rihoe + - riho 1 (46) 

where neo and nho are the homogeneous concentration of electrons and holes, respectively. 
This leads directly to: 

dp f _VE.\ 

^ = -^(;neoe-+n,oe . (47) 
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Together with Eq. (H5|) this leads to: 

(n.„e» + n„e-») {aP + PP^] = ^. (48) 
Substituting Eq. fH6|) into Eq. (!38|) . the Poisson equation for classical gas reads: 

dP 

— = -q l^rieoefeT - nhoe - + nhoj (49) 
This is a quadratic equation with respect to e^T. It has only one positive solution: 



—2^, ' 

and together with Eq. ( 148|1 it leads to the following equation: 

, ^'^/^--^ (51) 



Such equation in the context of the problem of polarization screening was obtained by Guro 
et al.— In the case of linear screening, the term ^ in Eq. fl5T]) can be neglected in comparison 
with rieo + n^o and Eq. ( 15T|1 can be rewritten in the form: 



The exact solution to this equation reads: 



-Po tanh ( |- ) (53) 



where 



g^(nco + nho)|«| 

In the case of a doped material, the concentration of minor carriers can be neglected, for 
example in the case of donor doping, Eq. ( I5^ can be rewritten as: 



5ci = W^^ (55) 

This spatial scale is about the Debye screening length for a linear media with a permittivity 
about 6f ^ 2TT/\a\. If the concentration of electrons inside the domain wall is much larger 
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than the homogeneous electron and hole concentrations inside the domain, Eq. ( 15T|) can be 
rewritten in the form: 

q oP/ox 

Changing variables P = Pqp, x = S'^ly, where: 

= 4^- (57) 
Equation fl56|) can be presented in a dimensionless form: 

The same characteristic scale S^^l was obtained by Krapivin and Chenskii^^. It corresponds to 
the nonlinear Debye screening length in linear media with permittivity about 27r/|a|, which 
is the typical length of charge screening in the case of classical statistics when homogeneous 
electron concentration is small in comparison with electron concentration in the screening 
area. The exact solution to this equation with boundary conditions p{—oo) = 1, p{oo) = —1, 
can be written as: 

The solution is presented in Fig. HI It can be approximated with the function p = 
— (2/71) arctan4?/. The function p = — tanh2?/ is given for comparison. 

This solution describes properly only the part of domain wall, where concentration is 
much larger than the homogeneous concentration. Thus the polarization profile given by 
Eq. f l59|) applicable to the case where the screening regime is nonlinear in most part of the 
domain wall. 



D. Screening with degenerate electron gas 

The density of states in valence and conduction bands for degenerate gas with parabolic 
spectrum can be obtained from the effective electron and hole masses as foUows^^: 
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FIG. 4. Exact polarization profile Eq. (j59p (solid line), approximation p = — tanh2y (dashed line) 
and approximation p = —{2 /it) arctan4y (dotted line). 



where rUe is the effective electron mass and mh is the effective hole mass. Below we neglect 
the difference in the densities of state in the valence and the conduction bands and between 
electron and hole effective masses. For numerical calculations we use the free electron mass 
m. 

For a degenerate gas, the Fermi functions can be approximated by a step function and 
Dirac-Fermi integral in Eq. (|5]) can be presented as 

FiMz) = ^z^l' (62) 

Two cases are possible in this approximation. Either valence band is fully occupied and there 
are electrons in conduction band, or there are holes in the valence band and the conduction 
band is empty. In the first case, using Eqs. f l60l) . fl62|) . Eq. (|5]) can be rewritten as: 

/(Mg^+Ep-Ec^ \ 

P = -g^ (63) 

In the second case, using Eqs. (16T1) . ( I62l) . Eq. ([5]) the charge density can be found as 
follows: 

P = \ -HoJ. (64) 
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Both Eqs. (163|) and (IMI) lead to the following result for dp/ dip: 

\ 1/3 

+ ^0 • (65) 



dp 3mq^ f \p\ ^ 



(37r2)2/3;j2 ^ q 

where hq is the homogeneous concentration of the majority carriers. 

Together with fH5l) and fl38l) this leads to an equation for the polarization inside the 
domain wall in the following form: 

1/3 o2 



3mg2 




dP 


(37^2)2/3/^2 




dx 



\ d^ P 

noj i^P + PP")=^- (66) 



In the case of linear screening, when the concentration of electrons is just slightly changed 
in comparison with their homogeneous concentration in the conduction band or the concen- 
tration of holes is slightly changed in comparison with homogeneous holes concentration in 
the valence band, Eq. ( 166|) can be rewritten in the form: 

{?,Ti^fi^n^ d^p 



^P + PP'- \ ; 1/3 g^ = o (67) 

This equation can be solved analytically, with the boundary conditions P(— oo) = Pq, 
P(oo) = — Pq, like in the case of uncharged wall, and the polarization can be found as 
follows: 



-Potanhl-^) (68) 



where 



3^ (-) 

Smq^^riQ \a\ 

is the typical half-width of a domain wall. This scale is the Thomas-Fermi screening length 
in linear media where the permittivity e is about 27r/|a|: 



^TF=oW .,,3 , 1/3 - (70) 

In the case of nonlinear screening, where the electron concentration in the conduction 
band inside the wall is much larger then the homogeneous concentration, Eq. (166!) can be 
rewritten in the form: 



'X 



2 



Changing variables: 



P = Pop, (72) 
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X = St^y, (73) 

where 



1/5 



^deg = ( 5 "Ji no ) , (74) 



Eq. (1711) can be transformed to a dimensionless form: 

3 1 d'p/dy' 

Thus, the characteristic scale for the domain wall half- width in this case is about (J^'g, which 
is of the order of the characteristic scale for nonlinear Thomas-Fermi screening for linear 
dielectrics with a permittivity of about 27r/|a|, i.e. the typical screening scale in the case 
where the electron gas in screening area is degenerate with the concentration much larger 
than the homogeneous concentration. The same scale (to within a factor of the order of 
unity) was obtained by Ivanchik~ for the case of internal screening in the monodomain 
state in a ferroelectric which is an intrinsic semiconductor. 

The exact solution of Eq. (1751) can be obtained in the following form: 

■"-\^) ioir^"^ ''''' 

The solution is presented in Fig. [5l It can be approximated with the function p = — tanhy. 

Similar to the case of nonlinear screening with classical electron gas, this obtained polar- 
ization profile applies to the case, where the screening regime is nonlinear in the most part 
of the domain wall. It is clear that the solution obtained above can be used also in the case 
of non-degenerate electron gas in the bulk of the domains, whereas inside the domain wall 
the gas is degenerate. 



E. Charged domain wall in a large sample. Applicability for different cases. 

Now, that the domain wall width for all regimes of screening is known, we can determine 
the conditions of applicability for each case. 

The electron concentration in the center of the " head-to-head" domain wall can be esti- 
mated as: 

p 

nc = + no. (77) 

where 6 is the domain wall half-width. 
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FIG. 5. Exact polarization profile Eq. ()76p (solid line) and approximation p = — tanhy (dashed 
line) . 



There are three main conditions, which subdivide the different regimes for domain wall 
screening. The condition of applicability of the formula for a strongly degenerate Fermi gas 
requires a small temperature in comparison to the Fermi temperature,— i.e. 



kT < ^— ^ 

2m 

In the opposite case the gas is classical. 

The condition for the linear regime of screening has the form: 



(78) 



n-nQ = Po/qS < uq, 



(79) 



hereafter riQ = max(ne,nh). 
If the obtained half-width 



5<5- 



neutral ; 



(80) 



the gradient term plays a decisive role and the domain wall half-width will be about the 

'^neutral '^'^c- 

The diagram presented in Fig. |6] was obtained with conditions Eqs. f l78|) -f l80|) . It is 
instructive to plot it in the coordinates of full- ionized dopants concentration instead of no. 
The zero doping corresponds to uq = 2ni, where rii is the intrinsic carriers concentration. 

For a non-degenerate gas, the line separating linear and nonlinear screening (lines 1, 2 in 
Fig. E]) can be described by the equation: 



no 



\a\Pl 
kT ■ 



(81) 
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FIG. 6. A diagram presenting areas that correspond to the different regimes of screening in a 
"head-to-head" domain- wall in a large ferroelectric crystal, depending on the ferroelectric energy 
density |a|-Po concentration of the fully ionized dopant. Here and N^, are the donor 

and acceptor dopant concentrations, respectively. Inside the region separated by dashed line, the 
domain-wall half-width corresponds to one of the scales obtained in Sect. IIIII Outside this region 
the gradient term becomes important and the domain-wall half-width is about the correlation 
radius re- 




For the degenerate gas, the regions corresponding to the linear and nonlinear screening 
are separated by the lines 3 and 4 in Fig. [HI and can be described by the equation: 

Lines 5, 6 and 7 in Fig. [6] separate the cases of screening by degenerate and classical gases 
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and can be described by equations 



37r2 
and 



1 /2kTmY^^ 

no = T^dcg = 1^ [—^ ) y^^) 



- ^.'-^s ^ (84) 

respectively. 

Lines 8, 9, and 10 separate the regions where the gradient term is important and can be 
described by equations : 

where Oq = h'^/mq'^ is the Bohr radius. 

The regimes of screening presented in this diagram work in extreme cases far from the 
borders of the corresponding regions. Near the borders, this regimes of screening will show 
a crossover behavior. Lines 2 and 4 on this diagram correspond to the full depletion. On 
these lines, the characteristic scales for linear and nonlinear screening can be matched, and 
the domain-wall half-width will be about the depletion width for the holes: 

5dep = — , (87) 

For the case of typical perovskites like BaTiOa and PbTiOs at room temperature, the 
domain wall corresponds to the regime of nonlinear screening by degenerate gas (see Sect. 

IMD. 



IV. FORMATION ENERGY OF A CHARGED DOMAIN WALL 

Now, that the polarization profile is known, we can calculate the energy per unit area of 
the isolated sample with charged domain wall, defined by Eq. fll3p . First we will neglect 
the energy associated with the surface boundary layer with atomic thickness. In this case 
the energy of the sample with charged domain wall are the same for "head-to- head" and 
"tail-to-tail" configurations. If the sample is large enough, the impact on the energy of the 
states with the domain-wall can be calculated separately for the domain-wall and the surface 
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adjacent regions of the sample for both configurations shown in Fig. [31 Then, the formation 
energy of the charged domain wall, Eq. f[T3|) . can be rewritten as follows: 

/OO j'OO 
[$h-h(a:)-$o(a;)]rfx+ / [^,^,{x) - ^o{x)]dx. (88) 
-OO J — OO 

where $h-h(3;) and $t-t(3;) are the free energies per unit area for the polarization profiles 
calculated with boundary conditions ( 15^ and respectively. 

Using Eq. ([T]), the free energy density, Eq. f[T^ . can be rewritten as: 

* = ^P^ + V + ^f^V + £iLp + ;,P»_ J^)V + t,,. (89) 



2 4 2 \dx J Stt V V '^a;^ 

We can neglect the electrostatic contribution in free energy in comparison with the ferro- 
electric one with respect to the small parameter eb/^^f- In our problem \P\ < \Pq\ that leads 
to \(3P^\ < \aP\, \K{d'^P/dx'^)\ < |KP/r^| = |2q;P|, and electrostatic term can be estimated 
as: 
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'aP + f3P' - K 



fd^p\ 


)■ 




— 


\ dx"^ J 




- 27r' ' 





and we can neglect it in comparison with aP^. Using 5 ^ rc we can also neglect gradient 
term | iii comparison with aP^. Finally, Eq. fl89|) can be rewritten as: 

$ = + ^-P^ + $eg. (91) 

For the case of a classical gas, the free energy density may be written as (see App. [B]): 
$(x) - $0 = (l^^' + ^P"" - ^i^o ) + (M^c - rih) - A;T(ne + - n,o - n^o)) • (92) 

Here, the first term corresponds to the lattice energy and the second one to the free kinetic 
energy of the electron gas. 

As shown in App. [Cl the energy of the charged wall linearly screened by classical electron 
gas can be obtained from Eqs. (IHHl) . (|92|) using solution Eq. (|53|) in the form: 

W,i = ^\a\P^S,,. (93) 



For the case of nonlinear screening by classical gas, using Eqs. fl56|) and fl92|) . the charged 
domain- wall energy can be presented in the form (see App. ICll : 



2 



W^cT = ^[E^ + In -8kT]], (94) 
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where is the density of states in the valence and the conduction bands. It can be shown 
(see App. [C]) that for Eg ^ kT, the second term in the parenthesis can be neglected in 
comparison with the first one, and the W^^ is given by: 

o p 

= -^Eg. (95) 

For the case of a degenerate gas, the free energy density can be presented in the following 
form (see Appendix [Bil : 



4 4 y V lOm 

(96) 

As shown in Appendix Oi the charged domain wall energy in the case of linear screening 
by degenerate gas can be expressed in the form: 

VTdeg = ^|a|Po'5deg, (97) 
and in the case of nonlinear screening 

IP 

W4 = 0.77|a|Po'< + -fE,. (98) 

It can be shown (see Appendix [C]) that for realistic values of spontaneous polarization 
and band gap (for i?g = 3 eV, Pq < 300 /iC/cm^), the first term in Eq. ( !98l) can be neglected 
in comparison with the second one, and the energy reads: 

9 P 

W^di = -^E.- (99) 

This relation is universal in the case of nonlinear screening for degenerate and classical gas. 
This is the energy needed in order to create sufficient electron-hole pairs for polarization 
screening. Although this result was strictly obtained for the case of a second order phase 
transition, one can show it is also valid for materials with a first order phase transition. One 
should note that if we neglect the surface energy the energies of configurations with charged 
domain wall are the same for "head-to- head" and "tail-to-tail" walls (see Fig. [3]). 

In the frame of the continuous theory the energy of the thin boundary layer can be 
taken into account introducing a phenomenological term 2^1 J^n in the formation energy of 
the domain wall, where is the polarization inside the domain, n is the normal to the 
surface directed toward the sample and Ci is the effective coefficient in the surface energy 
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(here we restricted our consideration with hnear term only). To obtain this coefficient the 
microscopical consideration is required. The difference between the energy of "head-to-head" 
and "tail-to-tail" domain wall is AdPofi. In the case of nonlinear screening the surface energy 
can be taken into account, introducing the effective changing of the bandgap: 



e: 



cff 



(100) 



for "head-to-head" and "tail-to-tail" walls respectively. The formation energy of charged 
domain wall with such notation have the form: 



W 



nl 



2Pr 



z-,cfr 



-e: 



g 



(101) 



V. SIZE EFFECT 



We will now address the question of the conditions for which the " head-to- head" config- 
uration is energetically favorable. Three other possible configurations: singledomain state, 
configuration with polarization equal to zero and multidomain state with antiparallel do- 
mains separated by neutral domain walls (see Fig. [7j) will be considered as competing sce- 
narios. Hereafter we will be considering the case of practical importance where the charged 
domain wall is in the nonlinear screening regime. Let us start from the competition between 

(a) . , , (b) 





(d) 
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FIG. 7. Schematics of polarization distribution for a state with charged wall (a), monodomain 
state (b), state with zero polarization (c), and the multidomain state with antiparallel domains 
separated by neutral domain walls (d). 



the "head-to-head" ("tail-to-tail") configuration and the paraelectric state by comparing the 
corresponding energies. We will take into account that for the sample with finite thickness L 
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the polarization can be different from the spontaneous one, due to depolarizing effect. While 
Eq. fllOip was obtained for the polarization equal the spontaneous one inside the domain, 
it can be modified for the case of arbitrary polarization inside the domain (if the screening 
regime in this case is also nonlinear). In this case the ferroelectric energy inside the domain 
should be taken into account. So that, for the sample much thicker than the domain wall 
width formation energy of charged domain wall can be estimated as: 

« , /5 d4 « D2^ , 2|Pd| ™ff 



W = L [^-Pl + ^P^ - -Po^J + , (102) 

where Pd is the polarization inside the domain. (Note that discussing the energy of a sample 
we apply the term energy to the difference between the energy of the sample and it energy 
of the material in the singledomain ferroelectric state times the volume of the sample.) It 
can be found as a polarization corresponding to the minimum of the energy Eq. (11021) . If 
E^^ > 0, the polarization in the sample with charged domain wall Pd < Pq- If Eg^ < 0, 
that is possible either for "tail-to-tail" or "head-to-head" configuration at high polarization 
effect of the surface, Pd > Pq. In a sample, thinner than some critical thickness we will 
find below, ferroelectricity can be completely suppressed due to the depolarizing field, and 
the configuration of zero polarization is energetically favorable. The energy (in the meaning 
specified above) of this way formed paraelectric state reads 

W = -L^P^ (103) 

Thus, from Eqs. ( |T3l) and f llOSp . the energy difference per unit area between the state with 
charged domain wall and paraelectric state can be written as: 

«r.2 , /5 , 2|Pd| ™fr 



AW = L \^-P^ + ^P^j + -^Ef. (104) 

The state with charged domain wall is energetically favorable if for some polarization Pd the 
difference AW < 0. This inequality has a solution if the sample is thicker than the critical 
value Lpar defined as: 

sVGEf 7.3 X Ef 
g|a|Po (l\o:\Po 

Thus, for sample thickness larger than Lpar the " head-to- head" configuration is more ener- 
getically favorable than the paraelectric state, otherwise the latter is favorable. 

It is worth noting such critical thickness is a/2 times larger than that that can be obtained 
from the condition 

Eth > Ef/{qL), (106) 
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where Eth is the thermodynamic coercive electric field. Such condition, written with the 
neglect of the surface poling effect was offered by Ivanchik— for the stability of a ferroelectric 
state with internal screening. 

Next we discuss the competition between the "head-to- head" ("tail-to-tail") configura- 
tion and the internally screened single-domain state. In the singledomain configuration the 
linear terms in surface energy is cancelled out and the energy of the internally screened sin- 
gledomain state equals a half of the energy of the configuration with charged wall, without 
surface energy term. In the case of nonlinear screening, one can find it as follows: 

W^on = y^g (107) 

(In writing this equation and further in the paper we cover only the situation in samples 
with the thickness L ^ Lpar where we can neglect difference between and Pq, thus we 
will use the energy given by Eq. (llOip .) This energy has the same meaning as Eq. ( !99l) . i.e. 
this is the energy that is necessary for the creation of sufficient electron-hole pairs for the 
polarization screening. Here an important remark should be made. In their treatment of 
the problem of monodomain state screening, Ivanchik^^ and Watanabe^ did not take into 
account the energy needed for the creation of electron-hole pairs. Therefore they obtained 
a much smaller energy for the internally screened monodomain state, leading to a much 
milder condition for the internally screened domains to be favorable in comparison with the 
domain screened by the charge in the electrodes. 

Comparing Eqs. f ll07p and f llOip one can find that one of the configurations with charged 
domain wall is energetically favorable if 2E^^ < Eg, i.e. > Eg. if the surface polarizing 
effect is not big enough, energy of the monodomain state is always smaller than the energy 
of the state with charged walls, thus the configuration with a charged wall is metastable, 
and charged walls might occur only as a result of growth kinetic of domains with elongated 
precursors (Fig. |8]). As a result of such growth either charged domain wall (Fig. |8]d) or 
lamella pattern with neutral domain walls (Fig. IHt) can be formed. Now let us discuss the 
competition between the "head-to-head" ("tail-to-tail") configuration and lamella pattern 
with neutral domain walls. In the case of lamella pattern with neutral domain walls the 
surface energy contributions are cancels out. The energy of the sample per unit area for this 
configuration is proportional to \/L and for both first and second order phase transitions it 
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FIG. 8. Schematics of possible domains growth processes. Domains growing in the form of elon- 
gated precursors (a), and two possible results of such growing: " head-to- head" domains (b) and 
multidomain state with neutral domain walls (c). 



W 477 X 0.26 X IKeutralV^' o fT nnQ^ 

VFd = 2 —= PoVL, (108) 

where and Ec are the lattice permittivity measured perpendicular and along the polar- 
ization, respectively, and IVncutrai is the energy per unit area of the neutral domain wall. 
Comparing energies given in Eqs. (11081) and (llOip . for E'^^ > one can find that the mul- 
tidomain state with neutral walls is more energetically favorable that the "head-to-head" 
configuration if L is smaller than: 

^-'i = (^] -A • (109) 

y q J 47r X 0.26 X Wneutral 

For a second order phase transition and near the transition temperature Tq, the neutral- 
wall surface-energy tends to zero according to the laws^ VTneutrai oc {Tq—TY^"^ and oc {Tq — 
T)~^, Ea is constant for uniaxial ferroelectrics and Ea oc (Tc — T)^^ for cubic ferroelectrics^. 
Thus tends to infinity and the multidomain state with neutral walls is always preferable. 
The situation is different for materials with first order phase transition, where close to the 
transition temperature, Lmd is finite. Here, depending on the sample thickness, either "head- 
to-head" domains or multidomain state with neutral walls is energetically preferable. Once 
one of these configurations is formed passing the phase transition, the spontaneous switching 
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between them, which is related with charge transfer, does not look very probable. 



VI. ISOLATED SAMPLE VS. SAMPLE WITH ELECTRODES 

The problem of the charged domain wall creation in the presence of metal electrodes was 
considered in a series of work by Ivanchik, Guro, Vul, and Kovtonyuk— and by Krapivin 
and Chenskiii^. After the mentioned authors, we will first neglect the surface energy effects, 
and will consider it separately. Below we show how some estimates for a sample with 
electrodes can be obtained, based on the results of the charged wall in the isolated sample. 
We will analyse the difference between these cases. 

We discuss the practically important case of nonlinear screening. In this case the energy 
of the charged domain wall is equal to the energy, necessary to create enough electrons and 
holes for screening. In the isolated sample, the only source of screening electrons is the 
valence band, and the creation of each electron- hole pair increases the free energy by Eg. In 
a sample with metal electrodes, a charge carriers exchange between the ferroelectric and the 
metal is possible, and the metal can serve as the source of the screening electrons or holes. 
The charge that will be created at the electrodes due to the electron transfer will screen the 
bound charge near the surface. 

To find the energy necessary for electron transfer we should take into consideration the 
chemical potential fi and the work function difference between the ferroelectric and the 
electrodes. For each electron that comes from the metal, the energy Eq — /i + — A{, where 
Ag and are the work functions for the electrode and ferroelectric respectively, should be 
added to the free energy of the system. The chemical potential for electrons is the Fermi 
level, and the energy of " head-to-head" domain wall can be estimated as: 

2P 

W = ^{Ec -Ep + A,- Af). (110) 
Q 

This relation can be clearly demonstrated with the band diagram from Fig. |9t^a),(b). 
When the Fermi level in metal is inside the band gap of ferroelectric (Fig. [ni(a)), the energy 
Eq. flllOp is positive, but it can be much smaller than in the case of isolated sample. In this 
case the charged domain wall is metastable, but with the energy smaller than in the case of 
isolated sample. If the fermi level in metal is above the bottom of the conduction band in 
ferroelectric (Fig. IHt^b)), the energy Eq. (IllOp is negative. The negative domain wall energy 
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means that the configuration with "head-to- head" wall is stable, while in the case of the 
isolated particle it is metastable. The same result with some additional small terms for the 
sample with electrodes was obtained by Ivanchik, Chenskii et al.— 
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FIG. 9. Band diagrams for different work-function differences between a ferroelectric and a metal 
when ferroelectric and metal are not in a physical contact. When the electrode and the ferroelectric 
are separated, the shift between the Fermi levels of the materials is defined by the work-function 
difference. The energy necessary for the creation of a unit free charge is shown with the bold 
arrow. This energy can be negative for " head-to- head" (b), or "tail-to-tail" configuration (d), and 
the domain wall is energetically preferable in this case. 



In the case of "tail-to-tail" wall, to provide screening, the electrons should move from the 
ferroelectric into the metal, the analog of the Eq. flllOp can be written as follows: 

2P 

W = —{EF-Ey + Af-Ae) (111) 
1 

When the Fermi level in metal is inside the band gap of ferroelectric (Fig. El^c)), the 
energy Eq. fill II) is positive, and the charged domain wall is metastable. If the fermi level 
in metal is below the top of the valence band in ferroelectric (Fig. El^d)), the energy Eq. 
(11 111) is negative, and the "tail-to-tail" domain wall is stable. 

The preferable domain- wall configuration ("head-to-head" or "tail-to-tail") is defined by 
the direction of the electrons transfer, which depends on the work- function difference between 
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the electrode and the ferroelectric. This result is quite different from the case of an isolated 
sample, where "head-to-head" and "tail-to-tail" configurations have the same energy. 

To take into account the surface energy term we should change Eg to E^^ in Eqs. f lllOp . 
f illip , and the result will depend not only on the electron affinity difference, but also on the 
surface energy coefficient. Depending on the sign of the electron affinity difference and the 
surface energy coefficient this effects can reinforce or work against each other. 

VII. NUMERICAL ESTIMATES AND COMPARISON WITH EXPERIMENTAL 
DATA 

Let us compare the results obtained with experimental data. First lets find what regime 
of screening corresponds to the typical perovskite materials at room temperature. For 
BaTiO^^ with Pq = 0.26 C/m^ ^ 7.8 x 10^ cgse and |a| = 3.6 x lO^^ the ferroelectric 
energy density |a|Po^ = 2.2 x 10^ erg/cm^ For PbTiO;?^ with Pq = 0.75 C/m^ ^ 2.2 x 10^ 
cgse and |a| = 0.05 one can calculate |«|Pq = 2.6 x 10^ erg/cm^. Both of this values are more 
than the critical value from Eq. (jHlj) for room temperature {2kTf^^m^/^/{3TT^h^) = 1.8x10^ 
erg/cm^. This means that the gas inside the domain wall is degenerate. Typically the homo- 
geneous carrier concentrations are non-degenerate. Thus, for typical perovskites, screening 
is provided by the degenerate gas in the nonlinear regime. The domain wall half-width can 
be found from Eq. f in|) and 6^1^ = 16 nm for BaTiOa and 5j'g = 2.7 nm for PbTiOa. The 
corresponding typical concentration in the domain wall ric ~ 10^'^ cm"'^ for BaTiOs and 
ric ~ 1.7 X 10^^ cm~^ for PbTiOs. These concentrations is still 2 — 3 orders smaller than the 
half-band-fill concentration and consistent with assumption we did to use the parabolic ap- 
proximation. At such concentrations the domain walls should have metallic conductivity^^. 
Interesting, no experimental data of enhanced conductivity in "head-to- head" domain wall 
is available. In contrast, in the recent experiment by Seidel et al.- the neutral domain wall 
show enhanced conductivity, while the charged domain wall in the same sample does not. 

We can also check, what conditions correspond to different regimes of screening. We use 
parameters of lead titanate^ with uq = 10^® cm~^ and temperature about 700 K (close to 
PbTiOs phase transition) to calculate the spontaneous polarization, corresponding to differ- 
ent regimes of screening. The dotted line on Fig. [TOT a) represents our material at different 
spontaneous polarization, corresponding to different temperatures. If the polarization is 
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small enough, i.e. the corresponding point is on the left from the point 1 in Fig. [TOT a). the 
regime of screening will be linear with classical electron gas. From Eq. flHTl) we can calculate 
that the spontaneous polarization corresponding to point 1 is about 8 jiC/cw?. Such values 
are possible close to the phase transition not in a bulk but in a thin film of PbTiOs, where 
transition expected to be of the second order—. From Eq. flM]) and Eq. ( 18 6 p one can find 
that point 2 in Fig JTOT a) corresponds to Pq = 30 ;uC/cm^ and point 3 to Pq = 2400 /iC/cm^. 
The last two values are independent of wq, we only assume that the homogenous concentra- 
tion is not degenerate. Point 3 can not be reached at any realistic spontaneous polarization, 
thus for typical materials charged domain wall is always wider than the neutral domain wall. 
This fact confirms assumption we did in Sect. IIIII The "head-to-head" domain wall with the 
thickness about the correlation radius may be possible in a materials with metallic carrier 
concentration, like GeTe^. 

Approaching phase transition we can pass through different regimes of screening, and this 
will be accompanying with the domain wall widening. The dependences of domain wall width 
from the temperature difference Tc — T are different for different regimes of screening and 
can be presented as 5 oc (Tc — T)~", where n = 0.7 for nonlinear screening with degenerate 
gas, n = 1.5 and n = 0.5 for nonlinear and linear screening with classical gas, respectively. 
Schematically this dependence is shown in Fig. [TOT b). 

We can evaluate the ratio between domain wall half-widths of charged and neutral domain 
walls using Eq. ( 174|) : 



(^neutral Pn /^fi:V2 U| VlO ' 



(112) 





where oq is the Bohr radius. For BaTiOs this leads to the ratio ^JJeg /"^neutral ~ 10. Equation 
(11121) is only slightly dependent on the parameters of the material, so that this ratio is 
virtually universal: charged domain wall width for typical perovskite ferroelectrics not too 
close to the phase transition should be about one order larger than the width of the neutral 
domain wall. This result is in good agreement with the experimental result of Jia et al.- 
where "head-to- head" and neutral domain walls widths were measured in a PZT thin films, 
using the negative spherical-aberration imaging technique. The width of a "head-to-head" 
domain wall is found to be about 10 unit cells, and width of neutral domain wall is about 
1 unit cell. As was already discussed, the results for the width of domain wall obtained in 
Hm is vahd not only for isolated samples, but also for a sample with electrodes, thus we can 
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FIG. 10. (a) A diagram presenting areas that correspond to the different regimes of screening in a 
charged domain-wah in a large ferroelectric crystal, depending on the ferroelectric energy density 
lajPg and the concentration of the fully ionized dopant. The dotted line corresponds to the fixed 
homogeneous concentration of carriers. Approaching the ferroelectric phase transition we can move 
from right to left along the dotted line, (b) Schematic of temperature dependence of domain wall 
width for temperatures below T^- 



compare it with the measurements done on thin films. 

We can also compare the energy of uncharged wall and the energy of charged wall, neglect- 
ing polarizing effect of the surface, in an isolated crystal at room temperature. The typical 
values of domain wall energies for different types of domain walls in BaTiOa and PbTiOs are 
presented in the Tab. HI The data for uncharged walls for BaTiOs is taken from Zhirnov^i, 
Hlinka and Marton^^, and Bulaevskii-^. For neutral domain wall of PbTiOs we used results 
of first-principles calculation by Foykko and Chadi^ and Meyer and Vanderbilt,— (at K 
the obtained neutral domain wall energy is about 100 erg/cm^) extrapolated to the room 
temperature using dependence VFneutrai oc (T — Tc)^/^ (see e.g.-^^.). To calculate the values 
for charged walls in BaTiOs we use Pq = 0.26 C/m^ ^ 7.8 x 10^ cgse and = 3.2 eV. For 
PbTiOs the spontaneous polarization Pq = 0.75 C/m^ ^ 2.2 x 10^ cgse and Eg = 3.4 eV. 

Now we can estimate minimal thickness at which the charged domain wall is energetically 
favorable in comparison with the paraelectric state. From Eq. fllOSp . without poling effect 
of the surface, one can find for room temperature for BaTiOa Lpar = 300 nm and for PbTiOs 
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TABLE I. Rough estimates for the domain wall formation energies BaTiOa and PbTiOs at room 
temperature. The estimates for the charge walls do not take into account the surface poling effect. 

Material Neutral Charged 

BaTiOs 10 1600 
PbTiOa 25 5000 



Lpsir = 80 nm. For the samples with thickness L < Lpar, zero polarization is favorable. Thus 
the existence of charged walls parallel to the film surfaces in thin films with the thickness 
about 10 nm (like e.g. in^) is possible only in the presence of injecting electrodes or strong 
polarizing effect of the surfaces. 

If the poling effect of the surface is small, the charged domain wall can be formed as 
a result of the domain growth kinetic at phase transition. The competing state, that can 
also occur as a result of such growing is the lamella structure with neutral domain walls (see 
Sect. |V]). To estimate minimal thickness, that multidomain lamella structure is energetically 
favorable in comparison with the charged domain wall at phase transition we should find 
energy of neutral domain wall near phase transition. It was shown that for second order 
phase transition VFneutrai oc Pq .— We will use the same dependence to set a rough estimate 
for the neutral domain wall energy for the case of first order phase transition .— The spon- 
taneous polarization at phase transition is about 0.16 C/m^ for BaTiOgr^ and 0.4 C/m^ 
for PbTiOa.^ Using domain wall energies at room temperature (see Tab. [T]) we estimate the 
neutral domain wall energies at phase transition PKieutrai ~ 2.5 erg/cm^ and PFneutrai ~ 4 
erg/cm^ respectively. Using Eq. fllOQp . neglecting poling ability of the surface, one can find 
for BaTiOg at the phase transition with Eg = 3.2 eV, ^ 5000, Ec ^ 2500^^ ^ 590 
/im. For PbTiOg at the phase transition with Eg = 3.4 eV, Ea ~ 1500, ^ 450^ Eq. ffT09|) 
leads to much smaler values L^^ = 80 /im. This means that for thickness from 80 /im to 500 
/im the charged domain wall is favorable in lead titanate, while in barium titanate laminar 
domain structure with neutral domain walls is. This estimate corroborates with the experi- 
mental results obtained by Surowiak et al.- Indeed, the structures formed during the phase 
transition in lead titanate and barium titanate crystals with thickness about 100 /im were 
studied. In the case of lead titanate "head-to- head" domain walls were observed, whereas 
in barium titanate an a-c domains laminar structure was observed where c-domains are 
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split into 180-degrees domains with neutral walls. One should also mention that the lower 
conductivity of barium titanate can be also a factor preventing the formation of "head-to- 
head" domain walls. Although in the case of nonlinear screening the homogeneous charge 
concentration does not affect the energy of the " head-to- head" configuration, the screening 
free carriers should be delivered to the domain wall and the higher conductivity can assist 
in the formation of the charged domain wall. 

In the case of the sample with electrodes, interesting observation could be done for ex- 
ample for BaTiOs with platinum electrodes.— >^ Neglecting poling effect of the surface for 
the case of the BaTiOs with platinum electrodes the formation energy of "tail-to tail" do- 
main wall is negative, i.e. domain wall is energetically favorable in this case. Indeed, 
neglecting poling effect of the surface Eq. f illip can be rewritten in a terms of band gap 
Eg and electron affinity x of ferroelectric in the form W = {2Pq/ q){Eg + x ~ ^c)- Using 
Eg = 3.2 eV, work function for platinum = 5.2 eV, x = 0-4 — 1-2 eV— >^ one can find 
Eg + X ~ = —(0.4 — 1.6) eV. This situation have not been observed experimentally. The 
possible explanation for this fact is the polarizing surface effect which works against the 
electron work-function difference. In the considered case of "tail-to-tail" wall this means the 
preferable direction of polarization due to the surface effect is out of the surface. 

VIII. CONCLUSIONS 

The full set of equations suitable for numerical calculations of a "head-to-head" domain 
wall in an isolated sample was considered. The effective boundary conditions for the case 
of the normal ferroelectrics (with e^/ei ^ 1) was introduced. Depending on spontaneous 
polarization and carriers concentration, different regimes of screening and different spatial 
scales of charged domain walls are singled out. For typical perovskites, not very close 
to the phase transition, this scale is about the nonlinear Thomas-Fermi screening length. 
For perovskite ferroelectrics this scale is about one order larger than the width of a neutral 
domain wall, which is in the good agreement with experimental results for PZT. Approaching 
phase transition by heating the sample we can pass through different regimes of screening; 
this will be accompanying with an appreciable widening of charged domain wall . The 
dependence of domain wall width from the difference between the temperature and the 
curie temperature are different for different regimes of screening. 
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TABLE II. Summary of charged domain walls characteristics for different regimes of screening in 
isolated sample (without polarizing surface effect) 



Regime Typical scale of domain Formation energy Remarks 

wall half- width per unit area of 

domain wall 

Half-width is about the Debye 
screening length for a linear media 
with Ef 2tx l\a\ 

Half-width is about the nonlinear 
Debye screening length for a linear 
media with Si ~ 27r/|a| 

Half-width is about the Thomas- 
Fermi screening length for a linear 
media with Si w 27r/|a| 

Suitable for the typical perovskites 
at room temperature. Half-width is 
about the nonlinear Thomas-Fermi 
screening length for a linear media 
with Ef ~ 27r/|a| 



Expressions for a domain wall formation energy at different regimes of screening was 
obtained. The information on the polarization profiles and domain wall energies without 
surface energy consideration for different possible screening regimes is summarized in Tab. 
Hn In the nonlinear regimes of screening this energy is equal to the energy necessary to create 
enough electron-hole pairs to fully screen the spontaneous polarization, being proportional 
to the electronic band gap of the ferroelectric. The results obtained are closely related to the 
problem of internal screening for monodomain state. It was shown that the energy of this 
screening is much higher than the energy obtained earlier by Ivanchik— and Watanabe^ due 
to neglect by the later an important contribution to it. In general, the surface contribution to 
the formation energy of charged wall can be comparable with the formation energy obtained 
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Linear bn\ = ^ ^ — — rr-r ||a|Pn<5ri 

classical 



Nonlinear 6^1 = '^E„ 

cl <7|a|-Po q 6 

classical 



Lmear 5,^^ = ^ _i_L_ -,\a\P,5^. 

degenerate 



Nonlinear = ( rX%j^ ) 

degenerate 



neglecting this effect. Either "head-to-head" or "tail-to-tail" domain wall can be favorable in 
comparison with the internally screened monodomain state if surface energy is big enough. 
If the charged domain wall is unfavorable, metastable " head-to- head" and "tail-to-tail" walls 
can be a result of the kinetic of domain growing. As a result of such growth either "head- 
to-head" domains or multidomain state can be obtained. Passing a second order phase 
transition, the multidomain state is always favorable, but for a first order phase transition, 
the favorable state depends on the thickness of the sample. 
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Appendix A: Boundary conditions. 



To check the self-consistency of the scale S, we should check that we can neglect first 
term in Eq. (12T|) . i.e. 



dE 

dx 



47r- 



dP 

dx 



using Eq. ([T]) this inequality can be rewritten as: 



ox OX'^ 



47r dP 

eb dx 



In solutions we obtained P < Pq, thus 3/3P^ < 3|a|P^ and using Sf ^ ey, we can get: 



{a + 3/3p2) 



dP 








StvOP 




Att dP 




< 








< 




dx 




ox 




Ef dx 


Eh dx 



To satisfy Eq. f lAip we should also show 

Q3p 



dx^ 



4n dP 

Eb dx 



For the case of linear screening, using solutions Pq = tanh(x/(5), one can get: 



KE-b d^P/dx^ 



KEb 



2sinh^(x/5) - 1 
1 + sinh^(x/5) 



Att dP/dx 27r52 
The maximum value of the function in parenthesis is 2, thus, using S > Ir^: 

KEbd^P/dx^ Eb /2rc^^ 



47r dP/dx Ei \ 6 
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-1^) <-«l. 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



In the case of nonlinear screening with the classical gas from solution (l59!l : 

l-p2)2(l-7/)|, (AT) 



An dP/dx tt5^ 

where —1 < p < 1. The maximum value of the function in modulus is 1, thus: 

< = — {^] < — < 1- (A8) 



Att dP/dx ~ n5'^ \ 6 J ~ e{ 

In the case of nonlinear screening with degenerate electron gas using solution ( I76|) . one can 
find: 

K£i,d^P/dx^ ^ K£b _ gb /rc\2 ^ £b . ..qx 

47r aP/ax 47r52 £f U 7 £f ^ ^ 

Appendix B: Free energy density for screening by classical and degenerate electron 
gas. 

In this section we will calculate the kinetic energy contribution in the free energy density 

$eg(n) - $egK)- 

In the case of the classical electron gas, from Eq. (I20l) for fc{p) ^ 1, the main term in 
the entropy density in the conduction band can be written in the following form: 

i 

where summation is done over all possible states in the conduction band. The Fermi function 
for a classical gas has the following form: 

fi = exp — , (B2) 



(B3) 



kT 

where Ei is the energy of the state. Together with Eq. ( IBip this leads to 
The energy density in the conduction band can be presented as follows: 

i 

Thus, the kinetic free energy density of classical electron gas in the conduction band can 
be found as follows: 

<^^^ = Eg - TSg = J2{EF + q^- kT) ff = (Ep + - kT) (B5) 
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On the same lines, the free energy density for classical electron gas in the valence band 
can be found in the form: 

^^^^ = -{E^ + qip + kT)n^ (B6) 

and the total kinetic free energy is equal to the sum of f IBSP and flB6P and can be presented 
in the following form: 

$eg = {Ey + q(p) (rie - rih) - kT{n^ + rih) (B7) 

Note, that Eq. (IB7p presents only the kinetic part of free energy, the electrostatic energy 
is already taken into account separately as ei,E'^/8TT in Eq. (ITS!) . The Fermi level is constant, 
thus in a sample with fixed number of electrons, Eq. (1B7I) can be rewritten as follows: 

$eg = q^^ine - nh) - fcT(ne + Uh) (B8) 

For degenerate Fermi gas, as was mentioned in section IIIIDt two situations are possible: 
either the valence band is fully occupied and there are electrons in the conduction band, 
or there are holes in the valence band and the conduction band is empty. With parabolic 
approximation the energy density can be calculated directly from Eq. ( IT8|) in terms of 
concentration, in a form which is universal for both situations as follows: 

^- = ^^l^<"='' + "'"> + "-^- 

The main term in the entropy density of a degenerate gas with respect to T/Tp, where 
Tp is the Fermi temperature, can be found from Eq. fl20|) as foUowa^S : 

^eg=(^)'''^fc^TnV3, (BIO) 

where n signifies either the concentration of electrons or the concentration of holes. 

The condition of applicability of the strongly degenerate Fermi gas approximation requires 
the temperature to be small in comparison to the Fermi temperature: 

kT < -n2/3. (Bll) 
m 

Thus, using Eq. (]B9[) 

TS^ « -^-^'^ (3) = 2 (B12) 
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and in Eq. (fT6l) we can neglect the term with entropy in comparison with the energy term. 
Thus, the free energy of the degenerate electron gas can be found from flB9p as follows: 

$eg = E,, = ^^J2'^\ nl/' + uT) + n^E,. (B13) 

So that for the degenerate gas we can finally get: 

$eg(n) - $eg(no) = 10 J (^e^' + ' " n'£) + (ne - n^o)E,. (B14) 



Appendix C: Calculation of the energy of a sample with "head-to-head" domain 
walls. 

1. Linear screening by classical electron gas. 

For a classical gas, free energy density is given by Eq. fl92|) . Using ( H6|) . the electron gas 
part of the free energy can be presented as a function of dP/dx: 

$cg(a;) - ^ego = qf ^"■eo - n^o - ~ (^eoe^ + ntoe'^^ (CI) 

where the potential (p can be found from Eq. ( l50i) 



kT 
Lf = — m 



IdP If IdP V . 

+ ^eO - ^hO + W + neO - ^ho + 4^60^^ 



2neo \ q dx V \ q dx 



(C2) 



In the case of linear screening, Eq. (1C1|) can be presented as an expansion with respect to a 
small parameter 9 = —, — and the second order terms should be taken into account: 

<?(neO+nho) ox 

$eg(x) - $ego = Axe + hrin^o + n^o)0^ + 0(6^) (C3) 

where A is some coefficient independent of 6, o{d'^) designates the higher order terms with 
respect to 9'^. 

Now we will obtain the relation between the ferroelectric and the electron gas parts of 
the free energy. With the change of variables: 

dP 

Eq. (152!) can be rewritten in the following form: 
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Integrating this equation once in P and taking into account that in the center of the domain 
t = and P = Pq, we can get: 

+ ~A T~ - 0-^7- — — ^ ■ y'^'^) 



2 4 4 2q^{neo + n^o) \dx 
Using Eqs. ( l92l) . ( ]C6I) and second order term from Eq. ( 1C3I) (first order terms will give zero 
after integration along the thickness of the sample with the boundary conditions, Eq. ( |37j) ). 
one can obtain: 

r°° f"" faP"^ BP^ aP^\ 
J [$h-hW-$o(x)]rfx = 2y [-^ + -^ f)dx. (C7) 

From Eq. (IC6p . |^ can be presented as follows: 



dP _ IfaP^ ^ f3P^ aP^ \ 2g2(n,o + nho) ^^^^ 



9x y V 2 4 4 y /tr 

With Eq. f ICSp . the integral in Eq. ( 1C7|) can be rewritten as an integral over the polarization 
as follows: 



(C9) 

In the case of linear screening by classical electron gas the polarization distribution for 
" head-to- head" and "tail-to-tail" walls are the same to within the sign, thus both terms in 
Eq. ( IHHi) are equal, and the total energy of a charged wall can be given as follows: 

W,i = ^\a\P^S,, (CIO) 



2. Nonlinear screening by classical electron gas. 

In the case of nonlinear screening, using Eq. f l55]) . Eq. (IBSp can be presented in the 
following form: 



$(x) - $0 = f + f - - V 





kT 


dP 


\dx J 


Q 


dx 



(cii) 



In the case of nonlinear screening by classical electron gas, Eq. (!56|) (with the variable 
change given by Eq. flC4p ) can be rewritten in the following form: 

«P + /3P^ = -— ^ (C12) 
q oP 
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Integrating this equation once, as we did for the case of the hnear screening, we will find 
the ferroelectric part of the free energy density in the following form: 



2 4 4 q \dx 

Then, the ferroelectric part of the energy of the sample, associated with the domain wall 
region can be found as follows: 

°° kT fdP\ ^ kT 2PokT l,_,^2eni 



From Eq. ( I46l) . in the case of the nonlinear screening, the potential can be presented in the 
following form: 



and the kinetic free energy of electron gas in Eq. (!92|) can be rewritten in the form: 

. r L (J-^) - 1) m = L (-?^) - 1) iP 

Q J-oo V V Q^eo ) ) \dx ) q J_p^ V V Q^co ) ) 

(C16) 

With dP/dx coming from Eq. (IClSp . Eq. (IC16P yields: 

Thus, the first term in (1551) can be obtained as a sum of Eqs. f lC14p and f lC17p as follows: 



-^fln» 1-4, (C18) 



q \ \ kTuco 

On the same lines, for the "tail-to-tail" domain wall region, one can obtain the following 
energy: 

The energy of the sample with charged wall can be obtained as a sum of Eqs. (IC18I) and 
f lUlQp . taking into account the mass action law = ngoJi-ho- 

,.1 «rF„ / f4\a\PS 



Using Hi = A^e-^s/^fcT^ ^j^^^g ^ VNcNy, Eq. (IC20|) can be rewritten in the following 
form: 



2Pa ^ , 4ATP„ / ^ fi\a\Pi 
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The second term in Eq. (lC2ip can be neglected in comparison with the first one, and the 
domain wall energy can be written as: 



2 P 

= -^E, (C22) 

Indeed, Eq. (]C21I) can be rewritten in the form: 

Wi = (^g + ^kT In - 8fcr) (C23) 

Typically in semiconductors is about several electron-volts and kT at room temperature 
is about 0.026 eV, thus 8fcT can be neglected in comparison with E^. Now we can estimate 
the term with logarithm in Eq. (IC2ip in comparison with E^. The density of states can 
be estimated as foUowa^S; 

The conditions for nonlinear screening by classical gas, corresponding to line 7 and to lines 
1 and 2 in Figj6]can be written as: 

FKo < p 



and 



^i<^o<^^, (C26) 
respectively. Using Eqs. ( 10251) . ( 1026^ . we can find following estimates: 

2kT In { -j-^ \ < 2kT ln(27r2) ^ QkT < E^. (C27) 

2fcTln (^^) » -2fcTln > -E,. (C28) 



Equations ( !C27ll and ( !C29ll imply: 



2A;Tln 



4|a|p2 



< Eg. (C29) 



kTN 

3. Linear screening by degenerate electron gas. 

In the case of a linear screening by degenerate electron gas, the amount of the electrons 
in the conduction bands remains constant and the carriers are just redistributed inside the 
band (either the valence band is fully occupied, and the electrons are redistributed into the 
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conduction band, or the conduction band is empty and the holes are redistributed into the 
valence band). Thus the electron gas energy can be found from Eq. flB14p as follows: 



This equation can be expanded with respect to the small parameter —jf- and the linear 
term can be neglected, because it will give zero after integration through the sample thickness 
due to the boundary conditions (1371) . Thus, Eq. ( ]C30I) can be presented in the form 



Integrating Eq. (1671) . one can calculate the local relation between the ferroelectric part of 
the energy and the part associated with the electron gas in the following form: 

In the linear regime of screening, the polarization distributions for the "head-to- head" and 
"tail-to-tail" domain walls are the same to within the sign, and with Eqs. (lB14p . flC3ip and 
flC32p domain wall energy can be found in the following form: 



oo 



With Eq. flC32p . the integral can be rewritten as an integral by polarization, and the 



domain-wall energy can be found as follows: 



4. Nonlinear screening by degenerate electron gas. 

In the most important case of nonlinear screening by degenerate gas, integrating Eq. ( ITTll 
the following relation can be obtained: 



2 4 4 hmq^l^ \dx J 

From Eqs. ( ]B13ll and ( 10351) . taking into account that J ricdx = 2PQ/q, the energy of the 
domain wall can be found as follows: 

<,.5rf^ + ^-^).. + ?^E. (C36) 



OO 



2 4 4 7 Q 
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where the polarization profile is obtained for the domain wall with the boundary conditions 
of Eq. fl39|) . As we did for the other type of walls, the integral in Eq. flC36P using Eq. (IC35P 
can be presented as an integral by P as: 

and Eq. ( 1C36I) can be rewritten in a following form: 

(r \ 2/5 pi p r^p 

- / (1 - z^I'dz X \a\PX., + —E, ^ 0.77 x \a\PXe, + —E,. (C38) 

The first term in Eq. flC40p can be neglected if 



Po «= ^El'' < {^y^' ^ lO^gse ^ 300/zC/cm2. (C39) 



V 1*^1 

Where we supposed Eg ^ 3 eV. Thus the energy of the domain wall can be written as 
follows: 

op 

W^lg = -^Eg (C40) 
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